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Abstract 

The goal of this paper is to discuss phase space analysis of some interacting Chaplygin gas models in 
General Relativity. Chaplygin gas is one of the fluids actively considered in modern cosmology due to 
the fact that it is a joint model of dark energy and dark matter. In this paper we have considered various 
forms of interaction term Q including linear and non linear sign changeable interactions. For each case 
late time attractors for the field equations are found. 


1 Introduction 


One of the active topics of modern cosmology it is an accelerated expansion of the large scale universe [I] - 
[5]. There are various approaches to solve the problem including dark energy [3], different modihcations of 
General Relativity m and other interesting ideas [TT] - [2] (to mention a few). There is an active research 
towards to all directions. If we keep General Relativity as a main theory of Gravity, then to describe 
dark energy we should assume either an explicit form of the EoS or the form of the energy density. It is not 
excluded, that the EoS of the dark energy can be a solution of an algebraic or a differential equation 
which could open a window towards fundamental description of darkness of the universe. One of the models 
actively studied in modern cosmology it is Chaplygin gas m, [iH], which has EoS of the following form 


Pc = Apc 


B 


( 1 ) 


where A, B and a are positive constants and pc it is the energy density of the gas. This is an example of 
the fluid which has an explicit EoS. It is well know that this fluid is a joint model of the dark energy and 
dark matter. Therefore, this is one of the main reasons to continue research on this model, despite to some 
critics m- Erom Eq. 0 the case of A = 0 recovers generalized Chaplygin gas EoS, and A — 0 together 
a = 1 recovers the original Chaplygin gas EoS. The best fitted parameters are found to be A = 0.085 and 
a = 1.724, while Constitution + CMB + BAO and Union + CMB + BAO results are A = 0.061 ± 0.079, 
a = 0.053 ± 0.089, and A = 0.110 ± 0.097, a = 0.089 ± 0.099 respectively [ID], [DT]. Other observational 
constraints on modified Chaplygin gas model using Markov Chain Monte Carlo approach found that A = 
0.00189l)J;i(5(^^^, a = 0.1079t)JJ^^^ at Icr level and A = 0.00189t!(;)!;(^)*^ with a = 0.1079l)J;29n at 2a 
level [22]. In recent Physical Literature different modihcations of EoS of Chaplygin gas could be found (see 
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for instance Hg - [17] and references therein). On the other hand, dark energy does give to an origin of a 
problem known as cosmological coincidence problem |28j . In this problem we should explain why 


^DM 

^DE 


= ro, 


( 2 ) 


where rg is a constant. To solve this problem an interaction between the components of the universe has 
been studied. An active research on interacting cosmological models reveals that interaction term Q can 
change the sing during the evolution of the universe [H], [30]. While before it was thought that it is not 
possible. In this paper we will consider cosmological models of interacting Chaplygin gas, where dark matter 
is a pressureless fluid. Our goal is to consider various forms of interaction term Q and to obtain late 
time attractors of the field equations describing the dynamics of such universe. Considered interactions are 
examples obtained from the forms presented below 


Q = q^{3Hbp + -fp), 


(3) 


and 

Q = 3Hbq^ , 


(4) 


where b, 7 , and n are positive constants, H it is the Hubble parameter, q it is the deceleration parameter, 
while p can be either the energy density of the effective fluid or the energy density of one of the components. 
Pi is the energy density of one of the components. We restrict our attention to the models where n = 0 and 
n = 1. Interactions described via Eq. (§ and Eq. for n = 0 have fixed sign during the whole evolution 
of the universe. Eor n = 1 interactions Eq. ([^ and Eq. Q represent sign changeable interactions, which is 
due to the deceleration parameter q. Interaction terms given via Eq. Q are ’’non linear” interactions with 
fixed (n = 0) and changeable (n = 1) signs, respectively. To study the dynamics of the large scale universe, 
in Physical Literature a very high accurate approximation is used, namely, we consider a mixture of the dark 
energy and dark matter as the content of the universe. Which does mean that the effective fluid is described 
as 

Peff = Pdm + Pde, (5) 


and 

Peff = Pdm + Pde- 


( 6 ) 


Therefore, with different assumptions about the dark energy and dark matter we are able to reconstruct the 
dynamics of the universe. Additional assumptions allowing us to construct various cosmological models we 
will present in section [^ 


The paper is organised as follows: In section [^ we will give the definition of the interacting models in 
modern cosmology and will present basics on the phase space analysis to find late time attractor solutions of 
the field equations of General Relativity. In section [^ phase space analysis is performed, late time attractors 
are found and classified according to their cosmological applicability. To save a place we presented only 
attractor solutions. Finally, discussion on obtained results are summarised in section]^ 


2 Interacting models and autonomous system 

It is well known that to describe the dynamics of the large scale flat FRW universe we need the following 
set of equations obtained from the field equations of General Relativity 


= 


SirGp 

3 


d 

a 


AtvG 

“iT 


(p + 3P). 


(7) 

( 8 ) 
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We suppose, that the cosmological constant A = 0, the gravitational constant G and c are constants with c = 
SttG = 1. On of the approaches to alleviate the cosmological coincidence problem is to consider interacting 
dark energy models. In modern cosmology an interaction between fluid components mathematically does 
mean [31] - [53] 

Pm + + Pm) = Qj (9) 


and 


Pa + ^H{pa + Pa) — —Q- 


( 10 ) 


The form of Q is determined under phenomenological assumptions, mainly, the dimensional analysis is used 
to construct interactions. It is reasonable to consider interactions which could improve previously known 
results and at the same time will not make the mathematical treatment of the problems complicated. It 
is widely believed that deeper understanding of the nature of dark energy and dark matter could give 
fundamental explanations of the phenomenological assumptions about interaction. It is known that a phase 
space of a dynamical system it is a space in which all possible states of the system are represented. To 
analyse the dynamical system of interacting Chaplygin gas we set |54] 


P<^ 

3772 ’ 

(11) 

Pc 


^ “ 3772 ’ 

(12) 

Pm 

^ “ 3772 ’ 

(13) 

= In a. 

(14) 


where a it is the scale factor. To have physically reasonable solutions we should have the following constraints 
0 < x < 1 and 0 < 2 ; < 1. At the same time we should remember that x and z according to Eq. Q should 
satisfy to the following constraint 

X + z = 1. (15) 

In terms of x and y EoS parameter of the Chaplygin gas reads as 


Wc = 


y 

5 

X 


while the EoS parameter of the effective fluid reads as 


^e// 


Pc 

Pc “f Pm 


(16) 


(17) 


because dark matter is considered as a pressureless fluid. It is not hard to show that the deceleration 
parameter q reads as 

9 = - l -= ^(1 + 3 ?/). ( 18 ) 

There is a huge number of articles presenting phase space analysis of different cosmological models |55j - 
[59j (to mention a few). Erom the next section we will start our study taking into account general algorithm 
of finding critical points of the autonomous system x' and y', where / it is the derivative with respect to 
N. Solutions of x' = 0 and y' = 0 should be found first, then the sign of the determinant and trace of the 
Jacobian matrix of x' and y' will point out the stability of the critical point. It is well known that if the 
trace of the Jacobian matrix is negative, while the determinate is positive, then the critical point is stable, 
because the real parts of eigenvalues are positive. This is in case of linear stability. From the other hand a 
stable critical point it is an attractor, which is we are looking for. Therefore, we need to find the range of 
the model parameters such, that to have a physically reasonable stable critical points i.e. 0 < Xc < 1 and 
0 < Zc < 1- 
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3 Phase space analysis 

In this section we will find late time attractors of the field equations for various interacting Chaplygin gas 
cosmological scenarios. For that we need to determine the form of the interaction term Q. Four different 
types of interaction are considered. Combining general experience known from Physical Literature we will 
impose the following constraints on the model parameters 


and 

In our calculations we used a fact that 


0 < a < 1 , 

(19) 

0 < 5< 1 , 

( 20 ) 

0 < 7 < 1 , 

( 21 ) 

A > 0 and B > 0. 

( 22 ) 

B 

^ 3m{Ax-y)' 

(23) 


3.1 Interaction Q = 3Hbp + 7 /? 

Interactions considered in this subsection are very well know and have been considered intensively. Our 
interest is to obtain the late time attractors of the field equations where the interactions are particular 
examples obtained from the following interaction 

Q = 3Hbp + 7 p, (24) 


where p could be either the energy density of the effective fluid or the energy density one of the components 


of the effective fluid. The general form of the interaction term given via Eq. (24) for two fluid universe reads 
as 

Q ^ 3Hb{pm + Pc) +^{pm + Pc), (25) 

which does allow us to obtain an explicit form of the autonomous system 


a;' = ^ = -3 (5 - 7 - J/(-l +x + y)), 


and 


y = 


dy 

dN 


= 3y{l + y)-3 


{-ay + A{1 + a)x) {b + x + y - 7(1 + y)) 


(26) 


(27) 


To find stable critical points (attractors), we need to solve x' = 0 and y' = 0 equations and determine the 
signs of the eigenvalues of the appropriate Jacobian matrix of x' and y'. Table does represent late time 
attractors for different forms of interactions obtained from Eq. (24). E.1.1 is a stable attractor for a, b, 7 and 


A satisfying to Eq.-s (19) - (22). EAA does represent Chaplygin gas dominating solution where ojc = — 1, 


We// = —1 and q = —1 i.e. ACDM model is recovered. Stable attractor EA.2 is a scaling solution since 


fie 


= 5 , 


(28) 


when 6/^0. Wc = —1 — 5 does mean that phantom phase is on face in this case. Attractor EAA is of the 
same nature as EAA, while attractor EAA has the same nature as EAA. Late time attractor E.1.5 has 
completely different nature, however in this case an accelerated expansion is not possible. E.1.6, E.1.7 and 
E.1.8 late time attractors are scaling solutions {b yf 0) since 


fic I — b 


(29) 
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Q 

S.P. 

X 

y 

SMbp^n ^Pm 

E.1.1 

1 

-1 

3Hbpc + 7 Pc 

E.1.2 

1 /( 1 + 6 ) 

-1 

3Hbpc + 'fPm 

E.1.3 

1 /( 1 + 6 ) 

-1 

3Hbpm + 'ypc 

E.l.A 

1 

-1 

SHbpm + 7 Pc 

E.1.5 

A—b—'y—r 

2A(1+7) 

6+7—^+r 

2(1+7) 

3Hbp + 7/9 

E.1.6 

1 - b 

-1 

3Hb{pc + Pm) 

E.1.7 

1 - b 

-1 

IP 

E.1.8 

1 - b 

-1 


Table 1: Critical points corresponding to various interaction terms obtained from general form Eq. ([^ for 
n = 0. r = yj 4A6(1 + 7 ) + (& — A + 7 ) 2 . 


In this case according to Eq.(18) 9 = — 1, while according to Eq. (16) 


Wc 


1 

1 - 6 ’ 


(30) 


EoS parameter of the effective fluid according to Eq. 0 is equal to —1. In Fig. 0 we have illustrated 
phase space portraits corresponding to E.1.1 and if.1.3. 




Figure 1: Phase space portraits for the models with Q = 3Hbpm + IPm and Q = iHbpc + jpm- 


3.2 Interaction Q = q{3Hbp + 7 /?) 

The second class of interactions that we will consider could be obtained from the following general form (two 
fluid case) 

Q = q{SHb{pc + Pm)+l{pc + Pc)) ■ (31) 
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Q 

S.P. 

X 

y 

q{3Hbpc + jpc) 

£ 2 . 2.1 

1/(1 - 6 ) 

-1 

iqiPc + Pm) 

£ 2 . 2.2 

1 

-1 

bq{3Hb{pc + Pm) + Pm+ Pc) 

£2.2.3 

1 + 6 

-1 

3Hbq{pm + Pc) 

£2.2.4 

1 + 6 

-1 

ZHbqpc 

£2.2.5 

1/(1 - 6 ) 

-1 

3Hbqpm 

£ 2 . 2.6 

1 

-1 


Table 2: Critical points corresponding to various interaction terms obtained from general form Eq. @ for 
n = 1 . 


It is know as sign changeable interaction, where the deceleration parameter q is used to have the desirable 


effect. For the general case given via Eq. (311, the critical points are solutions of the following two equations 


a;' = - (-6(1 + 3y)+j + y{-2 + 2a: + dq + Sqy)), 


(32) 


and 


y' = 32/(1 + 2^) + ^ (““2/ + "^(1 + (^(1 + 32/) + 2(a: + y) - {1 + y){l + 3y)-f). (33) 

Late time attractors for different interactions obtained from Eq. (31) are collected in Table Our analysis 
does show that late attractors are possible to obtain if 6 = 0 (Fig. ~§). 



Figure 2: Phase space portraits for the models with Q = q{3Hbpc + 'ypc) and Q = jqipc + Pm)- 


3.3 Interaction Q = 3Hb (^p + j 

In this subsection we will consider a non linear interaction. In Physical Literature there is an active discussion 
on non linear interactions (see for instance [60) . [61)1. The non linear interaction forms considered in this 
paper have a phenomenological origin and are based on the dimensional analysis. The general form of the 
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non linear interaction considered by us has the following form (for two fluid models) 


Q 


iHb 



Pij 

Pc “t" Pm 


( 34 ) 


where pij it is the product of pc and pm and if i j = m we have etc. In Table we summarize late 
time attractors which we have obtained using different forms of interaction term Q having similar forms with 
Eq. Q for n = 0. E.3.1 is a scaling late time attractor for 0 < a < 1, 0 < & < 1 and ^ > 0, because in this 
case 

^ = ^( 6-1 + n /1 + 6(6 + ^>)). 

It is not hard to see that Wg// = — 1, 9 = — 1 and 

_ _ 26 _ 

1 + 6 — A^/l + 6(6 + 6) 


(35) 

(36) 


Another scaling late time attractor is possible to obtain when 0 < 6 < 1/2, 0 < a < I and A > 0, if we 
consider the following interaction term Q 


Q 


iHb 



In Table it does correspond to E.3.2, for which 


rim _ 1 , Vl - 462 

~~2^ 2(1 - 2b)' 


It is not hard to see that 


- 26 

26 - 1 + \/l - 462 ■ 

Another scaling late time attractor is £1.3.3, when 0<a<l, 0<6<1 and A > 0 with 


(37) 


(38) 

(39) 


rtm 26 - 1 + v'l-4(6- 1)6 
OT “ 2(6- 1) ’ 

26 

"71-4(6- 1 ) 6 - 1 ’ 

and 9 = — 1, We// = — 1. Finally, if we consider interaction of the following form 

Q = 3i£6 (pm+ ) , 

\ Pc Pm ) 


(40) 

(41) 


(42) 


then late time attractor £1.3.4 could be found, which describes Chaplygin gas dominating phase with Wc = —1. 
In all cases an accelerated expansion is on face. Consideration of non linear interactions does give us 
qualitatively four different scaling solutions (among considered interactions). Therefore, future consideration 
of additional possibilities will enrich our experience. In Fig. ^ we have illustrated phase space portraits 
corresponding to £1.3.1 and £1.3.3. 


3.4 Interaction Q = 3Hbq (^p + 

Consideration of non linear sign changeable interactions is another phenomenological assumption of this 
work. Mainly, motivated by the results obtained in subsection |3.3[ we decided to consider the following 
possibility 

Q = 3Hbq ( Pc + Pm H - 7— 1 

V Pc + Pm / 
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Q 

S.P. 

X 

y 

3Hb {pc + pi/{pc + Pm)) 

E.3.1 

(-1 - b + ri)/2b 

-1 

^Hb i^Pc Pm “1“ Pm!^Pc “1“ PmT) 

E.3.2 

(—1 + 25 + r2) l^^b 

-1 

3Hb (pc + Pm + Pc/(Pc + Pm)) 

E.3.3 

(rs - l)/3b 

-1 

3Hb [pm + Pm! {Pc + Pm)) 

E.3A 

1 

-1 


Table 3: Critical points corresponding to various interaction terms obtained from general form Eq. Q for 
n = 0, where ri = \/l + 6(6 + b),r 2 = V1 — 45^ and = y/l — 4(6 — 1)6. 




Figure 3: Phase space portraits for the models with interactions 3i?& (pc + Pc/(/^c + Pm)) and 
3Hb (pc + Pm + Pc/iPc + Pm)) ■ 


We followed to the receipts of the previous subsections to find late time attractors for this case as well. 
Without going into all details of our calculations, we would like to report that we have considered two 
particular forms for the interaction term Q and found that late time attractors do not exist. Considered 
interaction terms are ^ 

Q = 3Hbq(p, + ^^), (44) 

\ Pc Pm J 

and ^ 

Q = 3Hbq (pm+ ) ■ (45) 

\ Pc T Pm ) 

We are not closing our study on this type of sign changeable interactions and we will come back to them in 
our next works. 

4 Discussion 

Subject of our interest has been a study of interacting Chaplygin gas cosmological models. Besides, we had 
a goal to construct some sign changeable non linear interactions. Instead to solve field equations for some 
initial conditions, we performed phase space analysis. Adopting appropriate notations allowed us to obtain 
an appropriate autonomous system and start a study on the stability of the critical points. Appropriate 




















restrictions from the cosmological and astrophysical studies know in Physical Literature are imposed on some 
parameters of the models to find physically reasonable solutions. Four types of interaction are considered. 
Two of them are fixed sing and sign changeable non linear interactions. For each case late time attractors 
are found. Our study showed that in case of some of the sign changeable non linear interactions, constructed 
in this paper, stable critical points are missing. While a study of the models when Q = q{3Hbp + 7/3) 
reveals that only one type late time attractor is possible, namely x = \ and y = —1 describing Chaplygin 
gas dominating phase with Wc = — 1 and <7 = — 1. The situation is completely different if we consider fixed 
sign interactions. In case of non linear interactions (among considered ones) four different types of late time 
scaling solutions were found, while for the models with Q = 3Hbp + 7/3 only two late time scaling solutions 
are found. One of them is described the universe in a phantom regime. More general cases of interacting 
models could be studied numerically and we hope to report some new results in our next articles. 
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